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In the preceding paper, the authors have presented an asymptotic solution for the 
field in the illuminated region of a large circular cylinder whose surface impedance around 
the periphery deviates from a constant value by a sinusoidal variation of small amplitude 
a. To 0(a), the reflected field comprises a specularly reflected ray and two first-order dif- 
fracted rays characteristic of a curved convex reflection grating. If the surface impedance 
varies ' 'slowly/' these three rays can be combined into a single specularly reflected ray 
having a reflection coefficient which depends solely on the local impedance at the reflection 
point. The "slowness" conditions necessary for the validity of this local reflection principle 
of geometrical optics are investigated and interpreted in physical terms. The results are 
presented in a manner which suggests their applicability to general, gently curved surfaces 
with slowly varying impedance properties. 



1. Introduction 

A plane electromagnetic wave incident on a large 
cylinder with a spatially periodic surface impedance 
around the periphery gives rise to a reflected field 
which may be interpreted as comprising a spectrum 
of rays appropriate to a curved convex reflection 
grating. These conclusions are presented in the 
preceding article wherein we have considered a 
circular cylinder whose surface impedance has a 
peripheral sinusoidal variation of small amplitude a 
superimposed upon a constant value [Marcinkowski 
and Felsen, 1962, henceforth referred to as III]. If 
the impedance varies slowly over a distance interval 
equal to the wavelength of the incident field, its 
periodic aspects in the vicinity of the specular 
reflection point lose their importance, and the 
specular and diffracted grating rays reaching a 
prescribed observation point are nearly parallel. 
From geometric-optical considerations it is rea- 
sonable to suppose that this ray bundle can be 
combined into a single specularly reflected ray 
having a reflection coefficient identical with that for 
a cylinder whose constant surface impedance is 
equal to the value of the variable surface impedance 
at the point of reflection. Such an assumption may 
be designated as the local reflection hypothesis of 
geometrical optics, and the conditions under which 
it obtains for the specular and the first order dif- 
fracted (grating) rays in III are discussed in this 
paper. These requirements are then phrased in a 
form which suggests their applicability also to im- 
pedance variations other than the one considered 
herein. 

To avoid unnecessary repetition, all suitable 
definitions are to be found in III. The reflected 
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rays considered here comprise a single, specularly 
reflected ray and two first order diffracted rays. 
Their analytical form is recalled in section 2, and 
the desired local reflection formula for a slowly 
varying surface impedance is presented. Section 3 
contains, in outline, a derivation of the conditions 
required for the validity of the local reflection hy- 
pothesis. Some physical interpretations of these 
requirements are provided in section 4. While these 
conditions^ are derived from the periodic impedance 
function Z(</>), the effects of periodicity in the neigh- 
borhood of the point of reflection play a minor role 
when Z(</>) varies very slowly. Hence, a formu- 
lation as in section 5, which expresses the require- 
ments for the validity of the local reflection 
hypothesis in a manner which makes no reference to 
the specific form of Z(0), may also be expected to 
apply to other slow impedance variations. 

2. Geometric Optical Fields to Order a 1 

To 0(a) in the perturbation solution, 2 the total 
field at a point in the illuminated region can be char- 
acterized in terms of an incident ray (the magnetic 
field vector H is assumed to be parallel to the cylinder 
axis), a specularly reflected ray of 0(a°), and two 
grating rays of 0(a). The axial magnetic field com- 
ponent is then given by 



67(p,0,<^ r )=6 r ^(p,^0 , ) + G r ^ (p,(/ > ,0 , ), 



(1) 



where the subscripts i and r identify the incident and 
reflected fields which are, to 0(a), 



<3 r < (p ) <k<£ / )=e~'* p0O8< *~*' ) 
6 r (pA,<i>')=G a T (p,4> ) 4>')+o{GKp,4>,4>','p) 

+QKp,<p,4>' ,-p)l 



(2a) 



(2b) 

2 a is the amplitude of the surface impedance variation, and all fields are ex- 
pressed as a power series in a. [see III, and Felsen and Marcinkowski, 1962, hence- 
forth referred to as I]. 
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Via the geometrical optics interpretation in III, 
the specularly reflected ray of 0(a°) is the familiar 
one associated with a constant impedance cylinder 



while a typical first order grating ray is 
G\(p,<t>,4>' ,p)=A Q {v)R,{p)D(p)< 



t iks(p) 



(3) 



(4) 



The A term gives the phase of the incident plane 
wave at the point of reflection, R and R x are the 
reflection coefficients, the D's are the divergence 
coefficients, and (iks) represents the phase along a 
reflected ray. We define the nonnegative integer 
p = 2ir8i/L, where a is the cylinder radius and L is 
the spatial period of the sinusoidally varying imped- 
ance. From this definition it follows that p/ka=\/L 
where k=2w/\ and X is the wavelength of the incident 
plane wave. All of these quantities have been pre- 
viously employed in III, eqs (6), (10), and (12). 

If the surface impedance of the cylinder is slowly 
varying, a consideration of figure 1, and of the grating 
law (7) of III, shows that almost every point in the 
lit region (0<^<7r/2) is reached by the bundle of 
three reflected rays comprising G T in (2b). In a 
typical situation, illustrated in figure 1, the three 
nearly parallel rays are reflected from the three dif- 
ferent points E, F, and G on the cylinder to the same 
point of observation P(p, <j>). Under the local reflec- 
tion hypothesis, it should be possible to combine 
these three rays into a single specularly reflected ray 
(q t = q x ) emanating from the point of reflection F in 
figure 1. This specular ray should have a reflection 
coefficient associated with a cylinder having a con- 
stant impedance equal to Z(9) defined by (2) in refer- 
ence III, where the angle $=#'±0*, Or, 5]0, is the 
value of (f> at the point of specular reflection, F. The 
value of Z(9) remains constant as the observation 
point P in figure 1 is moved along the ray FP, but 
changes for other locations of F on the cylinder. 
Under these conditions the local reflection coefficient 
RiiBi) has exactly the same form as R in (3) except 
for the substitution Z— >Z[l-f- a cos p(^>— </> )]. 
Therefore, we seek to combine the three reflected 
rays in (2b) into the single, specular ray, 



G T {p J <i>,<t> f )^A {i R l %)De :k \ 
having a reflection coefficient 



Ri(d i ) = 



cos Qi—Z\\-\-a cos p4>— O )] 



cos 6i-\-Z[l-\-a cos p(4>— </> )] 



(5) 



(6) 



Since (5) represents the result expected from con- 
ventional geometrical optics considerations applied 
to variable impedance surfaces [Keller, 1956], the 
conditions assuring its validity can be interpreted as 
requirements for conventional geometrical optics. 



3. Requirements for the Validity of the 
Local Reflection Principle 

Since a detailed derivation of the conditions for 
the validity of the local reflection hypothesis is 
somewhat lengthy, only an outline is presented in 
this section; the complete calculations are available 
elsewhere [Marcinkowski and Felsen, 1961, hence- 
forth referred to as II]. In the analysis, the geo- 
metric optical parameters for the grating rays are 
expressed in terms of small deviations from the 
corresponding parameters for the specular ray; for 
example, the path length GP in figure 1 is given by 
s(p)=s+As(p), with G falling below F for #>6. 
Thus, we define 

0<(p)=M-A0,Gp), 

O r (p)=e i +Ad r (p), 

s(p)=s+As(p), (7) 

and note that 6 r =di for the specularly reflected ray. 
From the geometry of the problem and the re- 
quirements of the grating law (7) in III, it may be 
shown that 



|A W | = 0(g)«l 



(8) 



uniformly in p and 0, provided that one imposes the 
restriction 



ka 



« cos0 f <l, O<0«<ir/2. 



(9) 



Therefore, these approximations break down near 
the shadow-lit boundary where di-^ir/2. (In this 
same transition region, the formula (2b) also be- 
comes invalid, see III.) Upon utilizing (7), (8), 
and (9), one may show that to 0(p/ka), 



Ae t (p)^T 
Ad r (p)^± 

As(p) c^— 



p 



ka cos 6i 2s + a cos d t 

p s-\-a cos Si 
ka cos B t 2s J r acosd i 

aAdj(p) sin 2 d t 



sindi+Adrtp) cos 6i 



=0 (10a) 



0,^0 (10b) 



(10c) 



where d r] =±e i if </>— <//^0. Since Adti+p) c~ 
— A6i(—p), it is evident that the ray bundle in figure 
1 is located symmetrically about the specularly re- 
flected ray (EFo^FG). For an observation point 
P(p, <j>) which approaches the surface of the cylinder 
(s->0), (10) shows that AS^p)^ and As(p)->0, 
while Ad r (p)-^O(p/ka)5*0. Consequently, even for 
the limiting condition of an observation point 
arbitrarily close to the surface, the bundle of three 
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Figure 1. Reflected rays to order a for a slowly varying 
sin face impedance. 

rays has a finite, nonvanishing angular deviation of 
0(p/ka). This result is a direct consequence of the 

grating law given by (7) in III. 

If the results of (10) are substituted into (6), (7), 
(12) and (23) in III, all the necessary grating ray 
parameters A , R u D, 6 i} r , and s may be evaluated. 
The calculations show that the fractional changes 
of the geometric optic parameters, in passing From 
the specular to the grating rays, are As(p)/s^Ar(p)/r 
^AD(p)/D^()(pfka). Use of these estimates and 
the results of (10) allows the transformation of (2b) 
into 

Griprf^'^AoRoDe^+autUR'De* 8 

r e ip(4>-<t>{)) e if(+p) JL e -ip(i-<t>o1 e if(-P) 1 (11) 

where 

J(p)c^kAs(p)+kaAdi(p) sin i ±pA6 i (p), 0,^0 

(12a) 



R'-- 



Z cos d t 



(Z+cosfl,) 2 



(12c) 



The restriction Re (Z)>0, already previously im- 
plied in the requirement for time-average power flow 
into the cylinder, has been imposed here to ensure 
the validity of the power series expansion of the de- 
nominator of Bi{p) hi (12) of III. In deriving (11), 
the amplitude of the perturbation parameter a has 
been explicitly related to the relative periodicity 
(p/ka) = (\/L) by the neglect of all terms of 0(a 2 ). 
This requirement implies that 



a0\ 



=0{a*) 



or 



ka 



(13) 



A similar estimate has been utilized by Hessel 1 1960] 
for an analogous study on an infinite plane surface 
with a sinusoidally varying surface impedance. If 



\m\<o(i)> 



(14) 



it is readily shown that (11) may be simplified to 
yield, to 0(a), the local reflection formula in (5). We 
conclude, therefore, that the requirement prescribed 
by (14) is a necessary condition for the validity of 
the local reflection hypothesis. 

4. Interpretation of Local Reflection 
Requirements 

The contributions to f(p) in (12a) arise from three 
differential effects which are not always of the same 
order of importance. The first two terms, of order 
unity, are much larger than the last term of 0(p/ka). 
The contribution described by the first term kAs(p) 
is due to the differential phase shift produced by the 
path difference As(p) of the different rays reaching 
the same observation point P(p, <j>) in figure 1. The 
second contribution, described by the kaA9 f (p) sin t 
term in (12a), arises from the A ) {j>) term which 
defines the phase of the incident plane wave at the 
point of reflection. Because of the cylindrical geom- 
etry, the phase of the incident plane wave at the 
three points of reflection E, F, and 6 in figure 1 is 
different, thereby introducing a differential phase 
shift into the reflected rays. The third contribution 
±pA6i(j)) in (12a) is due to the angle $(p) in the ex- 
ponential of (12) in III which provides the phase 
shift of the reflection coefficient R\(p). 

A consideration of the geometric optical solution 
given by (4) and shown in figure 1 indicates that the 
differential phase shift produced by the path differ- 
ence tends to cancel the differential phase shift pro- 
duced by the incident plane wave. This may be 
verified by making use of (10) in (12a). As a result 
of this cancellation, the first two terms of order unity 
in (12a) combine to form the first term of 0(p/ka) in 
(12b). In the very narrow back-scattering region 
where d^p/ka, this first term in (12b) is of O(0i) 
and therefore tends to vanish as 0<— >0. This agrees 
with physical expectations since the increments due 
to differences in the path lengths and in the phase of 
the incident plane wave both tend to be smaller in 
the vicinity of 6t = 0. (The apparent singularity pro- 
duced by the vanishing of the denominator does not 
arise since sin d r (p)^A6 r (p) as 6 t —>Q.) In this region 
the dominant contribution to j\p) is the phase shift 
term coining from the reflection coefficient R\(p). 

A 

This contribution depends upon the angle <t> at the 
point of reflection but is independent of the asym- 
metry (with respect to </>') produced by the phase 
angle <£ in the variable impedance Z{4>). 
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To obtain an overall picture of the behavior of 
j(p), it is desirable to make the following simple 
estimates, valid for most B i} 



Ad, 



cos 0*= 0(1) 

d r (p)=0(d i ) 



(15) 



In view of the approximate nature of these estimates, 
the requirement on f(p) prescribed by (14) seems 
overly precise and will be replaced by the simpler 
condition 



I/(2>)|«1. 



(16) 



Use of (15) in (12b) permits (16) to be written as 

L»2irZ , L >\/2ir (17) 



where we define 



L =2a\A6 i (p)\. 



(18) 



Z is the distance EFG in figure 1, measured along 
the surface of the cylinder, and can therefore be 
interpreted as the length of the cylindrical surface 
along which the variable impedance Z (<£) is sampled 
by the three rays reflected from the points E ) F, and 
6. The inequality (17) states that the local re- 
flection phenomenon obtains if the variable imped- 
ance interval L sampled by the diffracted rays is 
much smaller than the spatial period L. From (8) 
it is recalled that all the estimates in this section are 
subject to the basic restriction p/ka= X/Z <C 1. 
Therefore, if the sampling distance Z <X/27r, then 
(17) becomes ineffective and is replaced by the 
estimate Z»X. The relation (17) gives one 
equivalent of the restriction imposed by (16). 

Another and perhaps more useful statement may 
be obtained by substituting the relations for Ad t {p) 
and AB T (p)^ in (10a) and (10b), and the first two 
estimates in (15), into (12b). The inequality (16) 
is thereby transformed into a relation which ex- 
plicitly involves the wavelength X of the incident 
plane wave, 



Z»\(*s) 



1/2 



[a(3s+2a)] 1/2 
2s+a 



(19) 



This provides an estimation of the lower bound of the 
relative periodicity Z/X of the variable impedance in 
terms of the distance s=FP in figure 1 and the 
radius of curvature "a" of the cylinder. The im- 
portant behavior is given by the (ks) 1/2 term since the 
remaining terms are slowly varying. If (19) is 
evaluated for three different ranges of observation 
points s, one obtains the simpler relations 



Z»X(2£s) 1/2 



if — <s<Ca 
4t — 



(20a) 



/K \ 1/2 /K \ 1/2 

L^\(-ks) ~X(Jfan if s~a (20b) 

L^>\(jka) ' ifs»a. (20c) 

It is seen that i/X changes slowly from a depend- 
ence on the relative distance (ks) 112 in the near field 
s <C a to a dependence on the relative size of the 
cylinder (ka) 112 in the far field s> a. This behavior 
is in accord with physical expectations since only the 
relative distance ks would be expected to play a role 
near the surface, while the relative curvature (ka)~ l 
would be expected to become important far from the 
cylinder. If the point of observation is so close to 
the cylinder that s <C X/4T then the simple estimate 
L ^> X replaces the no longer effective estimate in 
(20a). 

The relation (19) involves explicitly the radius of 
curvature a. In the limit a^> oo } the c^ylindrical sur- 
face transforms into an infinite plane (19) then goes 
over exactly into (20a), thereby suggesting that 
this latter estimate should also predict the condition 
to be imposed on Z/X in order to have the local 
reflection hypothesis apply on a plane surface. For 
verification we note that (20a) agrees with a similar 
result obtained by Hessel [I960] who investigated an 
infinite plane with a sinusoidally varying surface 
impedance. This agreement highlights the negligible 
influence of the curvature of a large cylinder on the 
reflection phenomenon in the near field. 

However, the restrictions for the plane and cylinder 
differ significantly in the far field (s>>a). For 
the infinite plane, the length Z , along which the 
impedance is sampled by the three rays, becomes 
infinite as 5 becomes infinite. This follows from the 
geometry of the infinite plane and the requirements 
of the grating law. From the geometry it follows 
that the angle of incidence of each of the three 
incident rays is the same. From the grating law 
it follows that the angles of reflection differ by 
0(X/Z). Therefore, as we move a distance s along 
the specularly reflected ray, the sampling distance 
Z must continually increase so that, as s^>co, 
Lq^co for the infinite plane. 

The corresponding limiting behavior for the 
cylinder is quite different. As a result of the cylin- 
drical curvature, all of the angles of the different 
incident and reflected rays differ by 0(X/Z). As 
we move a distance s along a specularly reflected 
ray, it is only necessary that there be an angular 
deviation of the angles of incidence and reflection 
of 0(\/L) to change the direction of a grating ray 
by the amount necessary to satisfy the grating law 
and pass through the point P(p, 0) in figure 1. This 
angular deviation is easily secured by changing the 
point of reflection by the angular amount A6i(p) = 
0(\/L) which results in only a small fractional change 
on the circumference of the cylinder. Therefore, 
as s-^oo, the sampling distance for the cylinder 
L ^aO(\/L) which is a finite number as verified 
upon substituting the expression for Ad t (p) from 
(10a) into (18). This difference in limiting behavior 
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for the infinite plane and the cylinder is expressed 
mathematically by the different expressions (20a) 
and (20c). Evidently, if *>>a, the upper bound 
required by (20a) for the infinite plane is much 
larger than that required by (20c) for the cylinder. 
Therefore, far away from tin 4 surface, much more 
rapid changes in surface impedance are permissible 
for the cylinder than for the infinite plane, within 
the confines of the local reflection hypothesis. This 
is a direct result of the much smaller sampling 
distance L for the cylinder than for the infinite 
plane. 

It has been noted that as .9 varies from s<C<C.a to 
s>>a, the lower bound on L/X given by (20) 
changes from (ks) l/2 to (ka) 112 in a transition region 
where s^a. It is interesting to obtain the value s 
of s in this transition region in terms of parameters 
which do not involve the period L. Since (17) and 
(20) represent the same phenomenon, the lower 
bounds on L given by (17) and (20b) must be of the 
same order of magnitude. This provides the 
estimate 



s ^4t~ 



(21) 



The expression for s in terms of L and X has the 
same functional form as that for the distance from 
an aperture (having a diameter L ) to the region 
wherein the radiation pattern undergoes a transition 
from the (near-field) Fresnel to the (far-field) 
Fraunhofer type (a standard estimate in this case 
is *~2Z8/A). 

If the sampling distance L <C \/2w or if the length 
of the reflected ray s<X/47r then L and s can no 
longer provide significant estimates of the require- 
ments for the validity of the local reflection phenome- 
non. Under these conditions L and s are replaced 
by the wavelength X in the single estimate Z>>X 
wherein the wavelength itself provides the restriction 
on L because of (8) and (9). 

The rays reflected from the points E, F, and G in 
figure 1 also have a certain resemblance to rays 
propagating through a random, inhomogeneous me- 
dium. If a geometric optical approximation is to 
be valid for such propagating rays it would be 
reasonable to expect that there should be certain 
restrictions on the extent of the inhomogeneities. 
In making a comparison of these restrictions with 
our results, the distance of propagation in the 
medium is analogous to s while the range of the 
inhomogeneities is analogous to L. Chernov [1960] 
has considered the conditions under which a geomet- 
ric optical approximation is valid in a medium with 
spatially random inhomogeneities. It is interesting 
to observe that Chernov gives the conditions X<<Z 
and ^J\s<^<^L as necessary conditions for the 
validity of such an approximation. The first re- 
quirement is precisely the basic approximation p/ka 
<<0 introduced in (9) and used in (14). The 
second requirement is given by (20a) if we set all 
numerical factors equal to unity. 



5. Extension to More General Types of 
Impedance Variation 

Although the investigation herein is concerned 
with a cylinder having a periodic surface impedance 
variation, the effect of periodicity is obscured when 
p/ka<^l. Suppose that p = l with ka »1. In this 
instance the surface impedance as expressed by (2) 
in III goes through only one spatial period around the 
periphery of the cylinder. Viewed locally, one 
detects in any limited region of the cylinder surface 
an impedance variation devoid of any aspects of 
periodicity. It is reasonable to suppose, therefore, 
that the preceding restrictions delimiting the validity 
of the local reflection principle of geometrical optics, 
if properly rephrased, might apply as well to surface 
impedances with arbitrary, slow variations on a 
gently curved surface. The desired condition for 
local reflection should be phrased as a limitation on 
the allowed rate of impedance variation. Hence, 
we first determine the maximum relative rate of 
change of impedance occurring for the sinusoidal 
impedance variation Z(<£), for which the geometrical 
optics approximation retains its validity subject to 
restriction (19). 

Use of (2) in III £ or Z (</>) and of the relation a = 
0(\/L) in (13) yields 



1 dz{x) \ 
k <lx 



-"(!)" 1 



<1, 



(22) 



z(x) 



where z(x) is the variable surface impedance and x 
=a<t> is the distance measured along the surface; the 
subscript m denotes the maximum value. We now 
seek to express (\/L) in a manner which makes no 
reference to the sinusoidal impedance variation from 
which it has been derived. One such formulation 
utilizes the inequality in (19) which, when substituted 
into (22), leads to 



dz(x) 
dx 



z(x) 



< 



X (2s +a) 2 



s a(3s+2a) 

\ 

4tt' 



(23) 



s> 



The expression on the left hand side of this estimate, 
representing the relative impedance variation in 
an interval of one wavelength, is compared with 
quantities on the right hand side which contain the 
geometrical distance s along the reflected ray and the 
radius of curvature a at the point of reflection. 
Hence, this formulation of the restriction required 
for the validity of the conventional optics approxi- 
mation is phrased in a manner independent of any 
specific functional variation of z(x) and might there- 
fore be expected to apply to an arbitrary, slowly 
changing surface impedance on a gently curved 
surface. The restriction s > \/4t follows from the same 
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restriction in (20a). For s<\/4tt (23) is superseded 
by (22) in which the replacement 0(\/£) 2 -><<l 
is made for reasons analogous to those discussed 
in the previous section. It should be emphasized 
that z(x) is assumed finite in this expression and that 
the inequality becomes meaningless as 2^0, oo . As 
before it is of interest to examine (23) for various 
ranges of s/a. If we omit all numerical factors of 
order unity, then (23) becomes 



dz(x) 



dx 



z(x) 



s 



X 
<C - 

a 



if — <s<^a 

47T~" 



if s > a 



(24a) 



(24b) 



Evidently the allowed impedance variation is more 
rapid for s<Ca than for $>a. The reason for this 
behavior is connected with the negligible influence of 
surface curvature in the near field and its appreciable 
effect in the far field. 

The limited applicability of the local reflection 
hypothesis of geometrical optics to variable imped- 
ance surfaces stems from the fact that the reflected 
field is influenced not only by the surface properties 
at, but also in the vicinity of, the specular reflection 



point. For periodic variations, the diffracted rays 
provide a means of "sampling" the nature of the 
surface in the neighborhood of the point of reflection. 
For a surface with monotonic impedance variation, 
it is perhaps suggestive to construct a quasi-periodic 
equivalent for the region under consideration and to 
use the associated diffracted rays to sample the 
surface properties near the specular reflection point. 
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